Paraxial Hamiltonian for photons in two-dimensional photonic crystal microstructures 
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New solid-state physics based approach is developed for analysis of the paraxial light propagation 
in two-dimensional (2D) photonic lattices of coupled dielectric waveguides or microcavities. In 
particular, using Maxwell's equations, a non-Hermitian Hamiltonian eigenproblem with respect to 
the spinor wave function of a photon is obtained for energy-dissipating photonic microstructures. 
The Hamiltonian is suitable for almost the entire subclass of 2D structures encompassing arrays 
of semiconductor microcavities and microstructured photonic crystal fibers, characterized by light 
propagating mostly normal to the periodic lattice plane. Methods of numerical solution are discussed 
and the formalism is applied to a square array of coupled semiconductor microcavities, revealing 
reach possibilities for tailoring photonic band structure both in the photon energy and photon 
lifetime energy broadening domains. In particular, a feasibility to open a double photonic crystal 
band gap simultaneously in the energy and lifetime energy broadening domains is demonstrated. 
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I. INTRODUCTION 

Photonic crystal structures offer unique possibilities 
for controlling light-matter interactionsii^ii^ii^ and tailor- 
ing light propagation^ by introducing lattice defects^ii 
or photonic crystals heterostructure barriers;^ Among 
proposed so far photonic structure configurations, two- 
dimensional (2D) lattices are highly attractive in techno- 
logical aspect and have a realistic potential of finding ap- 
plications in novel optoelectronic devices and integrated 
photonic circuits. 

A particular subclass of 2D structures consists of ar- 
rays of coupled optical waveguides or microcavities. In 
these structures, a photonic mode propagates mostly 
along the waveguiding direction or cavity axis [vertical 
direction in FiglT] (a)], such that only a small lateral 
component kj^ of wave vector k undergoes Bragg refiec- 
tions. Such paraxial photonic structures employ lattices 
of period significantly exceeding the optical wavelength. 
Photonic crystal fibers^ are an example of implementing 
this concept. Matrices of phase-coupled vertical cavity 
surface emitting lasers (VCSELs) are another example of 
such 2D structures. ^^'^^ These structures allow photonic 
wave functions behaviour to be directly examined in peri- 
odic and quasi-periodic lattice a^'^i^^ and offer a possibility 
to implement photonic crystal heterostructures made of 
several photonic crystal materials with dissimilar band 
gapsi^l^ In addition, they might be loaded with optical 
gain or loss, which are show n^^i^^i^^ to impact band edge 
energies at high-symmetry points of reciprocal lattice. 

Many of these concepts related to photonic crystal 
structures originate from analogy to semiconductors (see 
for example, Refs. [l9l andl20l in the case of author). How- 
ever, despite an apparent simplicity of peculiar 2D pho- 
tonic crystal microstructures presented above, a simple 



approach utilizing a standard Hamiltonian framework of 
the solid-state physics has not yet been reported so far 
for photonic crystals. 

In this paper, a single-photon non-Hermitian Hamil- 
tonian is obtained for 2D photonic lattices of parallel 
dielectric waveguides or microcavities. It is derived by 
introducing a paraxial gauge transformation, which con- 
verts Maxwell's equations into a non-Hermitian Hamil- 
tonian eigenproblem with respect to a biorthonormal set 
of spinor wave functions for photons. (See Eq. ([50]) . the 
main result of this paper.) Properties of this Hamiltonian 
in case of lattices with inversion symmetry {e.g., widely 
used square and triangular lattices) are discussed and 
a method of numerical solution based on biorthonormal 
plane wave expansion is outlined. The numerical solu- 
tions are obtained for the structures with square sym- 
metry of the lattice, predicting a peculiar interplay of 
the photon energy and photon lifetime energy broaden- 
ing bands. These results envisage new possibilities for 
application of dissipative photonic crystals heterostruc- 
tures benefiting from the features of 2A^-dimensional 
confinement of photonic envelop wave functions in A^- 
dimensional photonic lattices. 

For the first time, a truncated form of this Hamilto- 
nian (adapted for the case of VCSEL arrays) has been 
briefly introduced in Ref. ^2^, without providing a discus- 
sion about its validity or proving orthogonality of its so- 
lutions. However, the symmetry and polarization struc- 
ture of the main lasing modes predicted by the model 
has been confirmed in Ref. [U by experimental measure- 
ments in square arrays of VCSELs. (See also Refs. [2^1231 . ) 
A k-p approximation of this Hamiltonian restricted to 
the eight bands^^ has been solved analytically, yielding 
the dispersion characteristics of the low-order bands in 
square lattices;^ A parabolic approximation obtained for 



the optical loss dispersion allows a single-band effective 
mass Hamiltonian to be introduced for modeling proper- 
ties of the main lasing modes in VCSEL arrays. Its finite 
difference element representation with the numerical grid 
step being the array pitch has the form of Coupled Mode 
Theory equations, offering an efficient tool for calculat- 
ing photonic envelope functions in complex quasiperiodic 
lattices.^ This method, which is detailed in the later 
Ref. has been used to interpret the experimental 
measurements of confined optical envelop functions in 
VCSEL-based photonic crystal heterostructuresi^ii^ 

In Ref,^, this Hamiltonian is used to analyze the 
Coriolis-Zeeman effect for photons in periodic lattices of 
microcavities, envisaging possibility of interactions be- 
tween photons and gravitational field. The basic steps 
of analysis in Ref. 2Q> might be used as guidelines for fol- 
lowing the treatment of photonic lattices in the present 
paper. However, note that in Ref. |20i, the dissipative 
effects are not taken into account. 

The paper is organized as follows. In Sec. [Hi using an 
equivalent cavity-unfolded representation of microcavi- 
ties, the Hamiltonian for lattices of microcavities defined 
by mirror reflectivity patterning is obtained. In Sec. IIIIl 
the Hamiltonian is extended to the case of microcavi- 
ties or parallel dielectric waveguides defined by periodic 
variations of refractive index. Sec. IIVI brings together 
the non-Hermitian Hamiltonian terms obtained in pre- 
vious sections, discuss the validity of such Hamiltonian 
approximation, outlines the method of solution based on 
biorthonormal plane wave expansion and reports numer- 
ical results of band structure computations in square- 
lattice photonic crystal microstructures. 



II. REFLECTIVITY-PATTERNED 
MICROCAVITIES 

The Hamiltonian for lattices of microcavities with pe- 
riodic variations of mirror reflectivity is detailed here 
on example of mirror-patterned VCSEL arrays (Fig. [1}. 
In Sec. Ill Al a typical structure of VCSEL array with 
reflectivity-patterned distributed Bragg reflector (DBR) 
is introduced. In Sec. IIIBl the constitutive equations 
in cavity-unfolded representation are obtained. Sec. Ill CI 
reports the model Hamiltonian. 



A. VCSEL array photonic crystal 

Arrays of vertical cavity surface emitting lasers are an 
example of photonic crystal made of evanescently coupled 
microcavities. These particular photonic structures uti- 
lize a usual VCSEL wafer incorporating a one- wavelength 
cavity with an optical gain medium (e.g. quantum wells) 
sandwiched between two DBRs [Fig[T] (a)] . The A-cavity 
might be regarded as a defect in the thick A/2-periodic 
DBR stack. In a solitary VCSEL, it defines the lasing 
mode wavelength, which fits the cavity roundtrip self- 
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FIG. 1: VCSEL array photonic crystal (a) Schematic illustra- 
tion of the wafer structure and metal-patterned top DBR com- 
position, (b) Cavity-unfolded representation of VCSEL array 
with equivalent (absorbing) g-layers representing the effect of 
reflections at the top (T) and bottom (B) cavity interfaces, 
(c) Simplified model system consisting of a Fabry-Perot cav- 
ity with reflectivity-modulated mirror, (d) Cavity-unfolded 
representation of the model system in (c). (e) Brillouin zone 
of the cavity-unfolded photonic crystal. A is the lattice pitch, 
a is the width of square VCSEL pixel. 



repetition conditions. In VCSEL arrays, the A-cavity 
fixes the main (longitudinal) wave vector component 
[along the z-axis in Fig[T]Ja)]. 

The photonic lattice of VCSEL array is defined by 
refiectivity patterning, R{x,y), of the top DBR.— The 
crystal unit cell consists of a highly reflecting VCSEL 
pixel surrounded by a grid of lower reflectivity. In the 
bottom-emitting structures, this is accomplished by de- 
positing a periodic pattern of metallic overlays on the 
top DBR, e.g.^ using gold pixels and chromium grid for 
pattern definitioniii In the top-emitting structures, the 
pixel positions and the grid are defined by air open- 
ings in the metallic film {e.g., gold film) deposited on 
the top DBRi^ The pattern is characterized in terms 
of a reflectivity contrast and lattice cell flU factor FF 
(the area ratio of the pixel and of the lattice cell) . The 
resulting composite DBR is of high reflectivity at the 
pixel positions (i?2~99.9%-=- 99.99%) and of low reflectiv- 
ity contrast between the pixels and the grid {2SR^1%). 
Since the roundtrip optical gain is low and uniformly dis- 
tributed across the VCSEL structure, such shallow re- 
flectivity pattern suffices to define position of the lasing 
microcavities. 

For a typical InGaAs/GaAs/AlGaAs VCSEL array 
emitting at 960 nm wavelength, the lattice pitch A is 
about 5 ^m, such that only a small transversal compo- 
nent /cj_~7r/A of the propagation vector k {k±/k ^ 0.03) 
undergoes Bragg reflections in the plane of periodic lat- 
tice. The angular spectrum of the most important low- 
order cavity modes is thus of narrow width ^ ~ A;_l /fc and 
centered about the main propagation direction (the z- 
axis direction). In what follows, the cavity modes are an- 
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alyzed in the paraxial approximation, within an accuracy 
of second-order terms ^^~10~^. The impact of reflectiv- 
ity pattern R{x, y) is taken into account as a second-order 
perturbation. (For a typical VCSEL array, the amplitude 
reflectivity variations assume 5R < 0.5 • fO^'^.) 

The optical gain is uniformly distributed across the 
cavity^I and has no influence on the cavity mode struc- 
ture. Uniformly distributed optical losses, such as the 
material losses (~ 0.1% in VCSEL arrays) and the out- 
put couphng losses at the bottom DBR (~ 0.16%), do 
not influence a curvature of the dispersion curves for the 
cavity modes as well. Their additive contribution to the 
optical mode frequencies and losses is ignored in the sub- 
sequent analysis and can be easily taken into account by 
shifting the frequency and loss spectra of optical modes. 

An optical mode of VCSEL array can be regarded as a 
superposition of standing waves coupled via Bragg scat- 
tering effects at the patterned mirror. The partial am- 
plitudes of standing-wave harmonics are defined thus by 
the boundary conditions at the cavity mirrors and not as 
a result of propagation within the structure, like in usual 
photonic crystal materials exhibiting periodic variations 
of the dielectric constant. To apply the usual method of 
orthogonal plane wave (OPW) expansion, an equiva- 
lent cavity-unfolded representation is needed. 

Multiple reflections at the cavity mirrors effectively 
translate the entire cavity into a structure that is peri- 
odic along the cavity axis (z-axis). The unfolded PhC is 
three-dimensional (3D) and can be analyzed in terms of 
propagating plane waves. Optical modes of a VCSEL ar- 
ray are thus represented by electromagnetic Bloch waves 
propagating in the equivalent 3D photonic crystal. 

A typical VCSEL wafer incorporates a few tens of var- 
ious dielectric material layers. To take into account a 
detailed composition of the wafer structure, the cavity 
has to be unfolded using the outermost layer interfaces. 
Thus reflections at the metallic overlay pattern in the top 
DBR and at the wafer substrate in the bottom Bragg re- 
flector have to be used for unfolding the VCSEL cavity 
[Figda)]. 

At each reflection from cavity output coupling inter- 
faces, the field amplitude reduces because of the waves 
leaving the cavity. In unfolded photonic crystal, the im- 
pact of subsequent reflections is reproduced by a periodic 
stack of energy dissipating layers (g-layers), each repre- 
senting the effect of a single mirror reflection [Figlljb)]. 
The zero propagation length between the incident and 
reflected waves at the output coupling interfaces of the 
cavity assumes that the equivalent g-layers are infinitesi- 
mally thin. Introduction of such thin absorbing layers 
allows both the energy dissipation (within the layers) 
and the continuity of electromagnetic field components 
(at each layer interface) to be taken by the model into 
account. 

In between the g-layers, the unfolded photonic crystal 
has the same sequence of dielectric layers as encountered 
in the VCSEL microcavity. However, such a detailed 
description of wafer structure composition is needed for 



analyzing technical issues, like, e.g., optimizing a spec- 
tral overlap between the gain peak and the cavity modes 
in function of the active region temperaturei^ This de- 
scription leads to intricate numerical model demanding 
lengthy numerical simulations. 

This paper is focused on general features of the light 
propagation behaviour in photonic lattices. Therefore, a 
simplified model system is used^i"— [FiglTJc)]. It consists 
of the A-cavity of VCSEL array in which the DBRs are 
replaced with ideal mirrors of effective DBR reflectivities. 
In this way, a Fabry-Perot cavity is formed in which the 
reflectivity of the upper mirror is modulated in two di- 
rections parallel to its plane. Such Fabry-Perot resonator 
replicates the boundary conditions for the field in VC- 
SEL cavity, reproducing thus the effect of fast oscillations 
within the periodic DBR stacks. Introduction of this 
model system is consistent with the Kapitza method^^ 
to analyze oscillations in a dynamic system. In this par- 
ticular case, it consists in separating fast longitudinal 
oscillations of the electromagnetic field (along the cavity 
axis) and its slow lateral oscillations (in the x-y plane). 

The model system treated here consists thus of a 
Fabry-Perot cavity of the length L=\/n with patterned 
reflectivity R{x, y) of one mirror [the upper mirror in 
FigU (c)] . The second cavity mirror is of uniform reflec- 
tivity and does not impact the cavity mode structure. In 
what follows, it is assumed to be perfectly reflecting. The 
dielectric material in the cavity is of uniform refractive in- 
dex n and impedance Z . The equivalent cavity-unfolded 
representation is a simple 2L-periodic (along the z axis) 
structure shown in Fig[TJd). 

Note that the periods (along the z axis) of the cavity- 
unfolded VCSEL array structure [Figlljb)] and its sim- 
plified model [Fig[ljd)] are different. As shown below, 
the absorbing g-layers have the same impedance as the 
neighboring dielectric materials such that Bragg reflec- 
tions in the z-axis direction cannot occur at the g-layers. 
Therefore, the period of the cavity-unfolded structure has 
no impact on photonic bands. 

Furthermore, the longitudinal wave vector component 
of a mode is set by the self-repetition condition at the 
cavity roundtrip. However, in the cavity-unfolded repre- 
sentation, it is not defined until such additional condition 
is superimposed, yielding kz = 2tt / L. 



B. Material equations for a cavity- unfolded 
structure 

There are many different conventions for the coordi- 
nate systems used to consider reflections, and the phases 
of reflection coefhcients are dependent on the coordinate 
system. Here, the following convention is used for the 
complex vector amplitudes of the incident and reflected 
waves: 
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where a is the reflection operator, Rs,p are the ampH- 
tude reflection coefficients, indexes s and p stand for po- 
larization perpendicular to the plane of incidence and in 
the plane of incidence, respectively. The convention ([T]) 
yields equal (in modulus and phase) reflection coefficients 
Rs and Rp at normal incidence, Rs = Rp{— R). Further- 
more, in the range of propagation angles S,=k±/kz consid- 
ered here, the angular dispersion of the DBR reflectivity 
is negligible, such that the approximation Rs=Rp is valid 
at oblique incidence as well. Indeed, for ^=0.03, the re- 
flectivity of a typical DBR decreases on Ji?~10~^. (The 
reflectivity amplitude at normal incidence is R^l.) 

The notion of absorbing g-layers in the cavity- unfolded 
representation of a structure is linked to the duality of 
symmetry operations.'^'' In the active notation, the re- 
flection operator tr in ([1]) assumes transformation of the 
vectors 'Ei(x,y,z) and H(a;,y,z). In the passive nota- 
tion this symmetry operation is applied to the coordinate 
axes. The effect of Ra on the flcld can be then regarded 
as a result of propagation through a layer with equivalent 
transmittance T=R. To match the phases of reflected (at 
the cavity mirror) and transmitted (by the equivalent g- 
layer) waves, the layer thickness has to be infinitesimally 
small. To impact the amplitude and phase of the trans- 
mitted waves in this limit, the material parameters of 
the layer should vary as oc with h being the layer 
thickness. 

In the cavity-unfolded representation, electromagnetic 
waves propagate in the periodic stack of absorbing g- 
layers [Fig[ljd)]. Material equations for this structure 
are obtained here by considering the amplitude of a plane 
electromagnetic wave in Fabry-Perot (FP) cavity with 
mirrors of uniform reflectivity. The cavity mirror located 
at z—0 is of reflectivity |7?|<1. As in Sec. Ill Al the sec- 
ond cavity mirror (at z=—L) is a perfect reflector. In 
the cavity-unfolded representation, the amplitude of the 
wave traveling towards the positive z-axis direction re- 
duces by a factor of R at each y-layer: 



E 
H 



Eq 

Ho 



(2) 



where to is the frequency of electromagnetic wave, Eq and 
Hq are the amplitudes of electric and magnetic fields at 
the coordinate origin z^O, 6{z) is the unit step func- 
tion e{z)^{0 (z<0), 1/2 (z=0), 1 (z>0)}, and j is an 
integer number (j'=0, ±1, . . .) enumerating g-layers (at 
z=0, ±2L, . . .) associated with the subsequent reflections 
at the cavity mirror. 

The electromagnetic wave ([2|) is written in the form of 
the Bloch wave composed of plane wave envelope function 
and periodic crystal part (the second exponential term is 
of the period 2L). The effective propagation vector k of 
the wave dH) is complex: 



Ini? 



(3) 



where um/c and r are the wave number and unit vec- 
tor along the wave propagation direction as seen from 



the cavity, and z is the z-axis unit vector. Eqs.([2l)-(l3]) 
show the evolution of wave amplitude and phase with 
successive roundtrips in the cavity. The imaginary part 
of k accounts for the decay of electromagnetic field in the 
cavity due to output coupling loss. The real part of k al- 
lows the cavity mode frequencies to be evaluated from a 
usual condition Re (2fc^L) =27rq with q being an integer. 

In the cavity [FiglTJc)], and hence in between the g- 
layers of unfolded photonic crystal [Figmd)], the electric 
and magnetic field components (O assume the relation- 
ships 



rxE = Zn, 



TxH = -Z-^E. 



(4) 



A correspondence between the electromagnetic energy 
flux in the cavity and in the equivalent unfolded struc- 
ture can be established assuming that the impedance of 
g-layer material is the same as in the dielectric material 
of the cavity. Otherwise, an impedance mismatch and 
Bragg scattering effects at periodic stack of g-layers will 
result in a retro-reflected wave propagating in the back- 
ward z-axis direction, having no counterpart in the FP 
cavity. The relationships (|4]) thus hold through the en- 
tire equivalent photonic crystal, including the absorbing 
g-layers as well. 

Substituting ©-Q in Maxwell' equations 



divB ^ 0, 
divB = 0, 



rofE 



rotH. 



1 dB 
1 dT> 



(5) 



within the accuracy of the ^^-order terms, one obtains 
the constitutive relationships for the periodic structure 
representing cavity-unfolded Fabry-Perot resonator 



D = eE + H X g 
..clni? 



g = -zz- 



UJ 



B = + g X E, (6) 
Y,5{z-2jL), (7) 



where e—n/Z and fi—nZ are the dielectric constant and 
magnetic permeability in the cavity, J^^S{^)d^=9{z). 
Eqs. ©-(IT]) assume the impedance matching through the 
entire structure, which thus exhibits no artificial photonic 
band gaps in the z-axis direction. 

Eqs. (in])-© are similar to the constitutive equations in 
a nonpermanent gravitational field induced by rotating 
coordinate framei^i^i^ Indeed, in Eq. ([T]), the reflection 
operator a at the cavity mirror includes the coordinate 
rotation ( a — /C2(z) where / represents the coordi- 
nate inversion and 6*2(2) represents rotation by tt about 
the cavity axis). The non-Galilean space-time metric— 
(g 0) within the layers represents the effect of mirror 
reffections and can thus be attributed to the coordinate 
rotations used to unfold the cavity. 

An apparent discontinuity of the electromagnetic held 
at g-layers [see Eqs. ([2]) and (jT])] is in fact a result of 
infinitesimally small g-layer thickness. Using Eq. ([7]) and 
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approximation 5(z)=lim/i^o 
one can obtain that 



, clni? _i 
-iz n 



(8) 



in a layer of small finite thickness h {e.g., located at the 
coordinate origin). It can be seen that the tangential 
fields E and H as well as the normal components of D 
and B are continuous at each interface of the g-layer, in 
agreement with the boundary conditions at moving (ro- 
tating) interfacesi^^i^i^Z, Within the layer, the tangential 
components of D and B vary as oc ft,^^, yielding a singu- 
larity in the limit /i— >0. 

The effective refractive index in the layer varies with 
the wave propagation direction^ 



20.33.38 



Tg- 



(9) 



For a plane wave propagating in the positive z-axis di- 
rection (z-x>0), the phase accrual in the layer reads 
/ yielding a complex number 



— (zg) dz = —i In i?, 

-h/2 c 



(10) 



where all other terms vanish in the limit /i— >0. At each 
g-layer, the wave amplitude thus reduces by a factor of 



0*0 



R, in agreement with Eq. 



The wave propagation in the positive 2-axis direc- 
tion of the structure ©-([T]) corresponds to multiple 
roundtrips in the FP cavity. However, this periodic struc- 
ture is nonreciprocal. Thus, the amplitude of a wave 
propagating in the backward direction (rz < 0) increases 
by a factor of (where |i?|<l) at each g-layer. This 
picture corresponds to an external electromagnetic field 
exciting oscillations in the FP cavity. The opposite z- 
axis directions in the cavity- unfolded structure ©-0 
are thus related by the time reversal operation. 

In the case of mirrors with uniform reflectivity, the vec- 
tor g in ([7]) shows the same transformation properties as 
vector g — irjxr composed of space-time components 
of metric tensor in the case of rotations. Namely, g is 
anti-invariant (g — g) under the coordinate inversion 
(P) and time reversal (T) operations. It can be seen 
that Eqs.([n])-(I7]) do not show any particular symmetry 
under the PT-transformation (coordinate inversion fol- 
lowed by time reversal operation). Therefore, the energy 
spectrum assumes complex eigenvalues, as opposed to the 
real spectrum in the case of pseudo-Hermitian systems 
that are invariant under the PT-transform.'^^'^*' How- 
ever, the PT transformation is important for the analysis 
presented here. Thus, in Sec. IIVI square-lattice struc- 
tures are treated using non-Hermitian Hamiltonian and 
biorthonormal set of wave function partners. The invari- 
ance of square lattice under (2D) coordinate inversion 
allows the partners of biorthonormal set to be defined as 
'^{t,x) and ^{t,x) = **(-t, -x). (In the case of dis- 
sipative effects, the results of substitution t —t and 



complex conjugation are different ,-4'- such that a wave 
function ^(^,2:) and its T-transform '^*{—t,x) do not 
coincide.) 



C. Hamiltonian 

The material equations ©-([Z]) define a cavity-unfolded 
periodic structure, which is equivalent to a Fabry-Perot 
resonator with mirrors of uniform reflectivity. The vector 
g in ([7]) does not vary with the position in the x-y plane 
(parallel to the cavity mirrors). Different from FP cavity, 
lattices of coupled microcavities, such as VCSEL arrays 
discussed here, utilize mirrors with reflectivity pattern- 
ing R{x,y). The cavity-unfolded representation of such 
mirror-patterned structures is assumed to obey the same 
relationships ([6]), however, with periodically varying (in 
the x-y plane) space-time coupling 



g 



, clnP(a;, y) 
-iz 



<5(z-2jL), 



(11) 



where the z-period of unfolded structure is 2L, and 
P(x, y) is the (amplitude) reflectivity of the cavity mir- 
ror [e.g., effective reflectivity of the top DBR in VC- 
SEL arrays). Note that in the model defined by Eqs.® 
and (|lip . like in the case of other photonic crystal struc- 
ture models, the field radiation effects at the boundary 
discontinuities^ are neglected. 

It can be seen that an electromagnetic Bloch wave 
propagating in photonic crystal with effective noniner- 
tiality pT|) is of the same form as indicated in Eq. ^ 
but the field amplitudes Eq and Hq vary periodically in 
the x-y plane. In the frequency w (real number) is 
the independent parameter of motion, while the Bloch 
vector k (complex) is defined from the dispersion equa- 
tion ^ and takes the effective propagation loss in the 
structure into account [last term in the right-hand side 
of Eq.®]. 

In fact, the independent variable parameterizing the 
dispersion curve can be chosen either as the wave num- 
ber or as the frequency of electromagnetic wave. In what 
follows, to obtain the standard Hamiltonian form of equa- 
tions, the vector k is used as a real- valued independent 
parameter. Respectively, the quantity hujq\^ is a quan- 
tized observable assuming complex values {q is the band 
index). The real part of observable ?iWgk is the photon 
energy and the imaginary part yields the lifetime broad- 
ening due to the optical loss. 

Maxwell' equations in a photonic crystal with periodi- 
cally varying noninertiality pip are solved here by sepa- 
rating the fast and slow field oscillations in, respectively, 
longitudinal (z-axis) and lateral (x-y plane) directions: 



Hqk 



^Ze,k(r±) 



(12) 



where the index q enumerates the photonic bands, the 
scalar function fqk{z) and the set of two vector functions 
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Gqk{r±) and hqk(rj_) are, respectively, the fast and slow 
oscillating Bloch wave components. In p^ . the field 
amplitudes e^k and h^k are normalized in vacuum. As 
shown below, such separation of variables is valid in con- 
ditions of the paraxial approximation and low contrast 
of photonic crystal lattice. 

As in the longitudinal part of the wave (|12p as- 
sumes a general form 



/9k (2) 



(13) 



where the first and the second terms are, respectively, the 
plane wave envelope function and periodic Bloch function 
part. The complex parameter c/igk defines the magnitude 
of periodic variations in the Bloch function amplitude 
and phase. This parameter is of the order of Ini? ~ 
Therefore, within the accuracy of the ^^-order terms, the 
wave (fT^ assumes the expansion: 



/gk(z) 



27r 



where 



2L 



(14) 



(15) 



The weak periodic modulation »7qk(-2) shows step-like 
variations at the positions of g-layers and thus repre- 
sents the effect of abrupt phase-amplitude variations at 
each reflection of the mirror. Note that the z-period av- 
erage values of the 2L-periodic function rjqk (odd func- 
tion) and its derivative drjqu^/dz (even function) are null 

(2x/-L^9k(z)rf2=(??gk)2L-0 ^ud {dijqk/ dz) :::iO) . By 
virtue of these properties, the longitudinal and lateral 
Bloch function components of electromagnetic wave 
can be analyzed separately. 

Substituting P^ - (|15p in Maxwell' equations for the 
curl of E and H, and taking the z-period average, one 
obtains the equations for the slow component propagat- 
ing in the lateral {x-y plane) direction: 



In R(x, v) , ^ „ 
^k. I 2 + 

In R(x, v) , ^ 

" 2L ' ^ + ^^ 



X egk=«-Wqkhgk, 



X hqk^-i-WgkGgk- 



(16) 



Then, by taking the difference between the z-averaged 
equation in (|16p and corresponding Maxwell' equation, 
one gets the equations for the fast oscillations in the z- 
axis direction 



drjqk 
dz 

dr]qk 
dz 



lnR^S{z -2jL) 
j 

InR^ 6{z-2jL)- 
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\nR 
~2L 



(17) 



z X h, 



gk = 



Eqs. p6)) can be converted into a Hamiltonian eigen- 
problem, provided that a photonic wave function 
■'/'«k(rj.) is introduced via gauge transformation of the 
fields e,k(i"_L) and \vq\^[vx_)- 

Wave states of a photon are usually expressed in terms 
of a three-component spinor wave function'^'^ (a photon 
spin is one). The spinor indexes represent polarization 
state, while the spatial distribution of spinor components 
accounts for the angular momentum. Thus, formally, 
three spinor components (with indexes 5^=0, ±1) have 
to be taken into account in Eqs. I|16p - p7|) . 

However, a state with 5^=0 cannot be realized in free 
spaced and only two spinor components (sz= ± 1) are 
independent. Thus, by virtue of the transversality of elec- 
tromagnetic wave in vacuum, only two field components 
might have independent spatial distributions. Obviously, 
the same holds for homogeneous dielectric media or pe- 
riodic photonic crystal structures. In the last case, the 
wave transversality condition is replaced by the coupling 
of the field components via Bragg scattering effects. 

Therefore, a photonic state is defined here in terms of 
two independent spinor components forming thus a two- 
component vector wave function xpq\s^(x,y). For parax- 
ial wave (|12p . it is convenient to choose the components 



i^qki^^y) ^^'^ V'gk (2;,y) of wave function ■ipq]s^{x,y) in 
the lateral {x-y plane) direction. A perturbation anal- 
ysis (not shown here) reveals that within the accuracy 
of ^^-order terms, the wave function components can be 
introduced in Eas. (|16p - (|17p via operator relationships 

Ggk = ^ • Ipqk, h,k = ^ • [Z X ■0gk] , (18) 

where (z-i/^gk) = 0, £ is the tensor operator 

iSa3 9 , 1 / 52 Sa0 



iSy)i 



kz dxp 2k1 \dxadx 



92 



K/3 



2 dx-ydx^ 



and twice repeated Greek indexes indicate summation 
over the x, y and z components. An explicit form of 
expressions (fT8 l) - ([T9| reads 



e(yk = 



A_L\ , 

l-ifcjjV'.k+.z 

/ A_L\ 

Kk=(l-^jzxt/;gk + 

V_L (V_L [zXt/^gk]) 

2fc? 



(V^L^/'gk) , V_L (Vj_l/J,k) 



2fc? 



., (V^ [zx-0gk]) 

?Z ; 



(20) 



Eqs. (|T8|) - (ll9p are used here as a definition of the gauge 
transformation introducing a wave function ■i/jgk- This 
gauge transformation is validated in Eqs. ((25)) and (|37p . 

In the gauge (|18p - p^ . the electromagnetic Bloch wave 
(fn)) reads 



qk 
rr(7) 
.■^9k. 



^ikzZ—iujt 



H-r?(z) 



27r 



4i\^,y), (21) 



where Cap-y is the completely antisymmetric unit tensor 
of Levi-Civita. The functions i/igk and zxipqy^ defines 
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the electric and magnetic fields of the main polariza- 
tion component. The first-order terms f-nd 

the second-order terms (oc-p-^^-^^^) contribute in lon- 
gitudinal and cross polarization components of the wave. 
Eq. (|2ip allows a complex polarization structure of inho- 
mogeneous wave to be taken by the model into account 
and is in agreement with the results obtained for Gaus- 
sian beams4i 

It must be noticed that ■0qk (and zxi/'gk) differs from 
the lateral {x-y plane) component of the field e^k (h^k)- 
However, the squared modulus ji/^gkCr^)!^ characterizes 
the energy flux in the z-axis direction. Thus, tak- 
ing a z-period average [Fig[Hd)] of the Poynting vector 
S=^Re (iE X H*)^^, one obtains Sz=-^ \i>qk\'^- In the 
paraxial approximation considered here, the energy flux 
in the z-axis direction is thus defined by the main po- 
larization component of the wave. The longitudinal and 
cross polarization components contribute to the energy 
flux in the lateral direction. 

As an example, consider a wave function i/^k = 
^^^ik^x+ikyy -^j^jj being a constant vector in the x-y 
plane. It defines a plane wave propagating in the direc- 
tion of k = (fca;, ky,kz). The energy flux associated with 
the wave is of the density Sz = -^ \ipo\^ along the z axis. 
The approximation (j2ip takes into account the transver- 
sality condition for electric and magnetic components of 
a plane wave. Thus, for the case of p-polarized (TM) 
wave (kj_ xi/'o=0), the electric field has nonzero compo- 
nents eoz=-k±/kz'ipo and eo_L=(l - jk^^/k^i/jQ, yield- 
ing the field amplitude eo—{l+^k'^/k'^)ipQ. The mag- 
netic field oscillates in the xy-planc with the amplitude 
h.o={l+jk^ / k'^)ipo ■ The energy flux in the direction of k 
is thus 5 = Soz{i+^kj_/k'^), in agreement with the vector 
calculus utilizing directional angle £,=k±/kz of the wave. 
By substituting e — > — h and h ^ e, a similar agreement 
can be readily proved for the case of s-polarized (TE) 
wave {^±1^0—0). 

Within the accuracy of the ^^-order terms, the inver- 
sion of Eq. ((20l) is straightforward, yielding the expres- 
sions for the wave function '0qk in terms of the fields e^k 
and z X h^k: 



^^^H'+Akih^-'^^ 2kr^' 



ZX'0,k = 



/ A_L\, .,(V_Lh,k) V_L(V_Lh5k 
= 1 1+— )hqk-«z ; — 



4fc 



(22) 



kz 



2fc? 



These relationships can be represented in the tensor op- 
erator form that reads 

where^ 630^637^=^^7 — SasS^^, and the inverse operator 
£ is defined by 



'a/3 



kz dxa 2k'^\dXadX[3 2 dx^dx^ 



(24) 



Thus, S is not the unitary operator and £ 7^ ^ • 

The paraxial gauge transformation (|2ip converts 
Maxwell' equations for the curl of E and H into the same 
form of a Hamiltonian eigenproblem with respect to the 



photonic state wave function "(/j^k = 1 (y) 

stituting expressions (|18p in pB]) and applying, respec- 
tively, the operators ezapSp}^ and f^J, in the first and 
second equations [in p^ ] one obtains 



cl%hiR{x,yy . , ^ , 
^ + ^ 1 ^gk = Tic^gkiAgk (25) 



2mo 



n 



2L 



where mQ=nhkz/c is the effective mass and p±=—iKS/± 
is the momentum operator in the lateral (xy-plane) direc- 
tion. The paraxial gauge transformation (I2ip represent- 
ing an arbitrary photonic state by a function xpq]s_{rj_) is 
thus validated by the fact that Maxwell' equations for E 
and H take the same form in this paraxial gauge. 

The first term in the Hamiltonian H [right-hand side 
of Eq. ((25)) ] is associated with the paraxial propagation 
along the cavity z axis and accounts for the dispersion of 
the longitudinal wave vector component in the dielectric 
material of the cavity {moc^ /n'^=hkzc/n). The in-plane 
kinetic energy (second term in H) and effective potential 
(third term) take into account the dispersion and Bragg 
scattering effects in the lateral direction due to periodic 
reflectivity pattern R(x, y) of the cavity mirror. The pa- 
rameter mo can thus be interpreted as the lateral effective 
mass of a photon in an empty lattice (in the case of R=\). 

The non-stationary Schrodinger equation for 
''/'<jk(i, r±) follows from by substitution id/dt lu. 
In the case of i?=l, within the accuracy of the time 
variable, the Schrodinger equation (PS)) is analytically 
similar to the scalar paraxial wave equation. How- 
ever, unlike a scalar field amplitude in paraxial wave 
equation, the spinor function i/j^k in Eq. (j25[) cannot 
be associated directly with any of the six components 
of electromagnetic field. The same remark applies to 
comparison between the Eg. ipS)) and the scalar 2D 
Helmholtz equation for microcavities."*^ 

The periodic potential U{x,y)=i^^ \nR{x,y) implies 
that i/?qk is a Bloch wave^^ composed of plane wave en- 
velope function and periodic crystal part 



■«/'gk(r_L) 



-Uqk(r_L), 



(26) 



where spinorial and angular parts of periodic function 
Uqk(r_L) allow the impact of photonic lattice symmetry 
to be taken into account in the analysis of electromag- 
netic field behaviour under the lattice rotations. Since 
the polarization anisotropy of mirror reflectivity R{x, y) 
is much smaller than the order of effects accounted for 
in the paraxial approximation [see Eq. ([T])] , the spin-orbit 
couphng term is neglected in Eg. ([25]) . All eigen states 
of the Hamiltonian ((25|) are thus doubly degenerate by 
polarization. As shown in Ref. IIO, this degeneracy can 
be removed by a symmetry breaking effects in nonperma- 
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nent gravitational field, when the photonic crystal rotates 
along the cavity z axis. 

Finally, the Hamiltonian H in Eg. ([25]) is non- 
Hcrmitian {H^^H). The general properties of non- 
Hermitian Hamiltonians were extensively studied in the 
past, yielding the conditions for a discrete real- valued 
spectrum of eigensolutions.'^^''^'' In Ref. 113, the interested 
reader can find a comparison between the eigenproblems 
of Hcrmitian and non-Hcrmitian Hamiltonians. In the 
case discussed here [Eg.ips])]. the most important results 
of these studies are related to the spectrum of H and 
orthogonality condition of its eigenfunctions. 

The Hamiltonian (|25p does not exhibit a pseudo- 
Hermiticity, as opposed to a class of PT-symmetric 
Hamiltonians, which are invariant under the time rever- 
sal followed by coordinate inversion. Therefore, the spec- 
trum of H assumes non-paired complex eigenvalues hujqk. 
Their complex conjugates can only be obtained in 

the spectrum of the adjoint operator . The real part 
of eigenvalues Iie{hLUqi^) is the photon energy and the 
imaginary part is the energy broadening due to a finite 
lifetime of photons in the cavity, Im(?iWgk) = h/2Tqk. 

The orthogonality of eigensolutions can be established 
using the biorthonormal set of functions^ consisting of 
the concomitant partners tpqi^ and xpqw of associated with 
the eigenproblems 

H \lPql,) ^hujqk \i>qk) (27) 

and 

i?+|t/i,'k')=K'k'l^9'k'), (28) 

respectively. Taking the difference between the matrix 
elements of Eg . ([77]) and the Hcrmitian adjoint of Eg. ipS)) . 
one obtains 

(na;,k-?.^,'k')(V^S'l^;kV0. (29) 

This relationship evidences^ the orthogonality of func- 
tions associated with different eigenvalues. In the case of 
Hamiltonian (|25[) . the spinorial structure of wave func- 
tions ipqii and xpqii as well as the plane wave envelopes 
and periodic Bloch function parts^ have to be taken by 
the orthogonality relationship into account, yielding 

(^S'I'/'S')= lirq,k,-Md'r^=S{k^-k'^)Sq,q. (30) 

Here, the integration runs over the entire x-y plane and 

(■^g'k'-V'<;k)='0g"k(2;,?/)*''/'gk^(a;,y) is the scalar product 
evaluated at a point r±={x,y). For photonic states at 
the same point of the 2D Brillouin zone (k^=k_L), the 
integral over the entire crystal ((30|) can be reduced to a 
lattice-cell integral^ 

(^SlV'S^)^^/ {rq-k-^,k)d'r^^Sq,q, (31) 



where f2j_ is the lattice cell area. For a state 
lipqk), the probability distribution function for the 
coordinates is thus given by the product i/'^k "'/'9k , 
as opposed to the usual expression |^gkP- In- 
deed the probability to find a photon at a point 

is (^£V(f-L-rl)|V£^)=iA;k(rl)-V',k(rl). In the 
Schrodinger picture, this expression allows the station- 
ary probability distribution to be associated with an 
eigen state. By contrast, the distribution I'^/'^kP varies 
in time as oc exp (— t/rgk)- This behaviour is expected 
by the photon hfetime considerations, since IV'gkP de- 
fines the the density of the energy flux along the cavity 
z-axis. Note that in the case of Hcrmitian Hamiltonian 
[Im(C/)=0], both distributions are stationary and indis- 
tinguishable. 

Substituting a Bloch wave ([26]) for the ■0gk and using 
the expression ■0gk = e''^^''^Uqk(r_L) for its concomitant 
partner, we obtain the orthogonality relationship for the 
periodic Bloch functions^ 

(4'^i"ik^) = ^ / («;'k-«,k)dv = v,- (32) 

The orthogonality of partners i/'gk and ipq'u' leads to 
the orthogonality relationship for the lateral components 
of the electric and magnetic fields (fT5|) . Since the gauge 
transformation operator £ [Eq. (|19p ] is independent of 
the material equations, the electromagnetic field associ- 
ated with the partner ■j/'gk is defined by the relationships 

eqk^S-ipqk, hqk^^'izxi/'gk]- (33) 

The scalar product in the integrand of Eq. ([50)1 can be 
then expressed in terms of the electric and magnetic fields 

Substituting Eq. for £ and using the relationship 

zx(v^(Vj.h;,k,)) = A^[zxh;,^,] - v^(v^[zxh;,k,]), 

one obtains that 

J {■4>*q'W^qi^)d^^±=- j z-[h;,k.xegk]d^r_L 

-eSV4[zxh;,k'])'"'})d^r^ (35) 
+ 2^/ (v^-jl^xh^k'KVi-e.k) 
-e,k(Vi-[zxh;,k'])})d'ri, 

where twice repeated index a indicates summation over 
the X, y and z vector field components. In the right-hand 
side of Eg. ([55)1 ■ the integrands in the second and third 
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terms are of the form (Vj^-F) with F being a vector 
field given by one of the expressions in curly brackets. 
As in Eqs. ([30|) - (|3T|) . the plane wave envelope function 
and periodic Bloch part-^ have to be accounted for in 
the electric and magnetic field components contributing 
to F. It follows that F can be represented as a prod- 
uct F=e'('^^^''^)'"^f (rj^) with f (rj^) being a cell-periodic 
function. 

Integration over the x-y plane can be then ac- 
complished by using the standard methods of the 
crystal field theory^S^ yielding /(Vj_-F)d^rj^ = 

S(kj_-k'j_)^ J^^^^{VjL-t)(fr±. Stokes' theorem allows 
the last integral over the lattice cell to be transformed 
into the contour integral along the cell boundaries yield- 
ing 

/ (V_L-f)rf^r_L= / V±x[zxi]-dS±=(f zx{-dl± 

Jccll JccU Jd{cc\l} 

where dS± and dl± are, respectively, the elements of the 
Wigner-Seitz cell and its boundary, and 9{cell} is the 
counterclockwise oriented contour in the x-y plane. 

The symmetry of the Wigner-Seitz cell assumes that 
for each point at the boundary, one can put in correspon- 
dence another point at the opposite boundary such that 
the two points are related by a primitive lattice transla- 
tion. The contributions from such points in the integral 
/a{ccii} 2 X f • dl± cancel out each other since at these 
points, the periodic function f (rj^) takes the same values, 
while the contour elements dl± are oriented in opposite 
directions. 

It follows that the second and third integrals in the 
right-hand side of Eq. p5|) vanish. Therefore, by virtue of 
the orthogonality relationship for the slowly varying 
components of the field reads 

y"z-[h;,k'Xe,k]rf2r^= -6ik^- k^) (36) 

where the integration runs over the x-y plane. Thus, 
within the accuracy of the paraxial approximation, the 
orthogonality relationship between the wave states from 
different bands [with indexes q and q' in Eq. pop ] is not 
influenced by the fast longitudinal (along z axis) oscil- 
lations of the electromagnetic field. Using p2)) . 
separating the integrals for longitudinal envelope func- 
tions and periodic parts [as in Eq. ([5^ ]. and noting that 
(^(?k)2L=0 with accuracy one obtains the orthogo- 
nality relationship for electromagnetic waves that reads 

J Z-[tll,^,XEgi,]d\^^^d{k,-k'J5{k^~k'^)Sg,g, 

(37) 

where integration is over the entire cavity-unfolded struc- 
ture. Note that relationship ([57)1 is in agreement with the 
waveguiding theories.— 

With the help of orthogonality relationship (PT|) . the 
solution of equations for the fast longitudinal oscillations 
of the field [Eq. (fT7|) ] is straightforward. As in the case of 



Eq. pS)) , it follows by applying the gauge transformation 

- —1 1 

and the operators —zx£ and 5 in the first and 

second equations in (fT7)) . Within the accuracy of 

order terms, both equations yield the same form of the 

master equation for the periodic part 77gk('Z) of the fast 

longitudinal wave component 

=^(l'gk2^Hz-2jL) - —, (38) 

i 

where '/>ijk~C^ (or less) is the diagonal matrix element 
(/)gk=— *(V'gk^l lii^lV'gk^)- The off-diagonal elements have 
smaller magnitudes and, within the accuracy their 
contribution cannot be accounted for in the right-hand 
side of Eq. p8)) . In this case, the left-hand side of Eq. p8)) 
is also null {{'4''g?l\^p'qZ^)dr]qk/dz=0 for q' ^ q). 

The integral of Eg. ([38)) is the expression given in (fT5|) 
but now with the parameter (f)q\^ expressed as the diag- 
onal matrix element of effective crystal potential. This 
verifies the separation of variables in the paraxial gauge 
transformation ((2T|) . 

The matrix elements (j)qk provide a way for more ac- 
curate formulation of the paraxial approximation con- 
ditions in the Hamiltonian ([^5)) . Thus, the effective 
crystal potential in ([^S]) has been assumed so far to be 

of the relative order that is f^^<C^^. Us- 
ing the roundtrip self-repetition condition in a A-cavity 
fcz=27r/L, the condition of low-contrast patterning can 
be represented in the form 

\lnR{x,y)\<A7re. (39) 

However, as evidenced by Eq. (|38p . it is more important 
that the effective potential variations are small in aver- 
age, at the scale of the crystal lattice cell. For a particu- 
lar photonic band, the last condition can be expressed by 
evaluating the diagonal matrix element in (j39p . yielding 

|(^;^)|lni?(x,^/)|^J^))|<4<^ (40) 

where the matrix element is the same as \(t>q]i\ in Eq. (j38p . 

The condition (j40p verifies the paraxial approximation 
(PS)) within a particular set of photonic bands. For struc- 
tures with simple lattice cell topology (like in Fig. [1]), 
this condition is less restrictive. For such structures, 
further simplification is possible in the most important 
case of low-order photonic bands, which are typically 
characterized by relatively smooth wave functions. In- 
troducing the reflectivity contrast parameter 6R as a 
measure of the maximum variations in lni?(a;, y), the in- 
traband matrix elements ([^0]) can be then estimated as 

\{ij'gl>\\iiR\Tpl^l')\^FF\6R\ for the lattices with smaU fill 

factor FF<^ {e.g., for the lattices with small-size pixels 
in the case of structures shown in Fig[T|) . In the oppo- 
site case of large fill factor FF>^ (large-sized pixels in 
FiglT]), the matrix elements are '^(l — FF)\dR\. 
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FIG. 2: Schematic illustration of the paraxial photonic crys- 
tal structures defined by periodic variations of refractive in- 
dex, (a) Arrays of coupled index-guided microcavities {e.g., 
etched VCSEL arrays^^) or parallel waveguides (like in the 
cladding region of photonic band gap fibers^^). (b) Comple- 
mentary, antiguided-array structures of microcavities defined 
by negative-index step (e.g., antiguided VCSEL arrays^) 
or low-index core antiguides (like in photonic liquid crystal 
fibers^). A is the lattice pitch, a is the pixel size. 



It follows that the paraxial Hamiltonian (PS)) is valid 
for the structures with reflectivity patterning contrast 



\5R\ < i'^An/FF 



< ^2 



iFF<l/2), 
(Ht:/{1~FF) (FF>1/2). 



(41) 



In most practical cases, optical microcavities satisfy this 
condition. 



III. DIELECTRIC LATTICES DEFINED BY 
REFRACTIVE INDEX VARIATIONS 

The approach developed in SecHIl is focused on lat- 
tices of optical resonators realized by mirror reflectivity 
patterning (Fig[T|). For completeness of the Hamiltonian 
(PS)) . it has to be extended to similar lattices of coupled 
microcavities or parallel waveguides defined by periodic 
variations of the dielectric constant (Figl2]) . Examples of 
such coupled microcavities with positive [Figl^l^a)] or neg- 
ative [Figl^ljb)] refractive index contrast are, respectively, 
etched VCSEL arrays^ or matrices of antiguided VC- 
SEL resonators fabricated using a two-step organometal- 
lic chemical vapor deposition (OMCVD) growth. As 
shown below, the approach based on Hamiltonian for- 
malism (j25p applies equally well to arrays of paral- 
lel waveguides^ [Figl^fa)] and low-index-core antiguides 
[FiglDJb)]. These photonic crystal materials are employed 
in the cladding of photonic band gap fibers^ and pho- 
tonic liquid crystal fibers 

In this section, the Hamiltonian (|25p is thus extended 
to almost entire subclass of 2D paraxial photonic crystal 
structures characterized by the light propagation mostly 
normal to the periodic crystal plane. As it can be ex- 
pected, the Hamiltonian for such dielectric lattices [see 
Eq. (|52p ] is analytically similar to the scalar paraxial wave 
equation but differs from it by the use of two-component 
spinor wave functions for photons. 

Note that a general form of the Hamiltonian 
for dielectric lattices can be obtained directly from 
Eq. (I25p . assuming that R=l and introducing the 



refractive index variations n(x,y)~^ — {n~^) and ef- 
fective mass mo=hkz/c{n~^) in the first term of 
Eq. ipSl) . These yield the effective periodic potential of 
U^moc'^{^){-^^-{^)). However, such calculus does 
not reveal the maximum refractive index contrast Sn sat- 
isfying the conditions of paraxial Hamiltonian approxi- 
mation. As shown below, the allowed refractive index 
variations assume the condition (|59[) , which for low-order 
photonic bands, can be represented in the form 

< ^/FF {FF < 1/2), ^^2^ 

<^/{l^FF) {FF> 1/2), 

where FF is the lattice cell fill factor, A is the lat- 
tice pitch, A is the wavelength (in vacuum) and ^ = 
(ri~^)A/2A is the order parameter (f<d) corresponding 
to a propagation angle of state at the boundary of the 
first (2D) Brillouin zone. Unlike the second-order varia- 
tions SR in reflectivity-patterned structures, the refrac- 
tive index variations dn can thus be treated by the model 
as first-order perturbations [compare Eqs. (|lTl) and (|^ ]. 
As discussed in Sec. IIV[ this allows the paraxial Hamil- 
tonian to be applied to some of the structures exhibiting 
high contrast of refractive index variations ((5n~l), as in 
the case of holey photonic crystal fibers^" or arrays of 
micropillarsi^ 

To analyze separately the impact of refractive index 
variations on photonic band structure, we assume that 
g=0 in constitutive equations (O. In the case of cou- 
pled microcavities, this suggests that the uniform loss 
and gain distributions are neglected and the reflectiv- 
ity of the output coupling mirror is R=l. As in the 
case of mirror-patterned structures, the vertical cavity 
composition is accounted for by using an effective (com- 
plex) refractive index profile^ n{x,y). The Bragg scat- 
tering effects set in, conditioned by the periodic pat- 
tern e(x,y)—n{x,yY / II. For the structures considered 
here, the magnetic permeability /i is assumed constant 
(/i = 1 in the optical spectrum range). Thus, in ac- 
cordance with the relationships e{x, y) = n{x, y)/Z{x, y) 
and fi = n{x, y)Z{x, y), both the refractive index n{x, y) 
and the impedance Z(x,y) exhibit periodic variations. 

The cavity-unfolded representation (Sec.|TT|) can be ap- 
plied to the matrices of dielectric microcavities as well. 
In this particular case, it effectively translates the cavity 
into a structure that exhibits a translational symmetry 
along cavity z-axis, allowing thus the correspondence be- 
tween an array of microcavities and equivalent structure 
of parallel dielectric waveguides. Electromagnetic waves 
propagating in such waveguiding structures (in the z-axis 
direction) can be expressed using the paraxial wave ap- 
proximation (jl2p with separated fast (longitudinal) and 
slow (lateral) wave oscillations. 

Since R—1, an electromagnetic wave propagating in 
a cavity-unfolded structure does not show sudden am- 
plitude variations at g-layers [see Sec. Ill B| . Therefore, 
as in the arrays of parallel waveguides as in the cavity- 
unfolded lattices of dielectric microcavities, the fast Ion- 
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gitudinal component of the wave is a uniform plane wave 
showing no periodic modulation in the z-axis direction 
[see Eq. (fT4l) ]: 



??,k(^) - 0. (43) 
The electric and magnetic field components then read 



Z(r^)-i/2h,k(ri) 



(44) 



After substitution of ([44|), Maxwell's equations for the 
curl of E and H yields the equations for the slowly vary- 
ing wave components in the lateral xy-plane direction: 



n{x,y) 
c 

n{x,y) 



ifc^z-|-Vj_- 
ifez-|-V_L- 



2Z{x,y) 
2Z{x,y) 



xh, 



9k 



(45) 



In order to convert these equations into the Hamiltonian 
eigenproblem, an appropriate gauge transformation has 
to be applied for the fields e^k and hgk- A general form of 
the gauge transformation, which takes into account the 
variations of Z{x,y) and n{x,y) = fi/Z{x,y) and, at the 
same time, encompasses the particular case of a uniform 
dielectric medium (pTt . reads: 

e,k = • ^,k, = Z-''£Z'' -zxip^k, (46) 

where, as in l|2ip . '0gk is a two-component spinor wave 
function, the tensor operator £ is defined in and 
the parameters a and b will be adjusted to fit the spinor 
transversality condition (z'0gk) — 0. 

The corresponding inverse transformation reads 

t/.,k=^"^"'^""-e,k, zx-0,k=^"''^"'^'-h,k, (47) 



where £ is the tensor operator introduced in Eq. 
Note that 1/n can be substituted for Z in Eqs. (|i5)) - ((i7)) 
since for considered here structures, the magnetic perme- 
ability fi is constant. 

In these expressions, the spatial derivatives of Z (or 
1/n) should match the paraxial approximation condi- 
tions, which thus limit the contrast of dielectric ma- 
terials composing the lattice. Introducing the average 
impedance (Z) and variations 6Z{x,y)=Z{x,y)~{Z) , one 
should assume that the relative variations dZ{x,y)/{Z) 
are of the order of ^. In that case, the gauge transforma- 
tion (P51) - (|T7l) allows spatial derivatives of SZ(x, y) to be 
taken into account by the model: 



5Z{x,y) 



< 



VSZix,y) 



k,{Z) 



<e 



(48) 



Otherwise (if SZ{x,y)/ {Z)^£^'^), in conditions of the 
paraxial approximation, the effect of periodic variations 
6Z{x, y) cancels out in the gauge transformation (146^ - 
(|47p . yielding thus expressions for a uniform dielectric 



material [Eqs. (fT8|) and ([23]) ]. The conditions in terms of 
refractive index variations Sn{x, y) follow from Eq. (j48p 
by substitution Z since Z{x,y)^fj,/n{x,y). 

Following the analogy between photonic crystals and 
semiconductors, it is interesting to observe that the ratio 
between two expressions in pS)) reads 



VSZ 



kJZ 



< 



\76n 



kz6r. 



< 1 



(49) 



This corresponds to the conditio n^°i^^i^° of smooth po- 
tential variations | ^JJ'^ | ^1 in the effective mass method 
widely used in solid-state physics. In the case of paraxial 
light propagation considered here, the condition (|49|) ver- 
ifies the use of effective mass mQ=hkz / c{n~^) in analysis 
of photonic bands in the lateral propagation direction. 

By substituting the gauge transformation ((46|) 
and applying the operators —ilie^apZ~^£'^^Z'' and 

ihZ°'£j^^Z^'^ in the first and second equations (|45p . both 
equations are converted into a similar form with respect 
to spinor functions zxi/'qk and rpqi^- 



cTikz 
n{x,y) 



zxi/jqk - {n 
1 \ , 



chAj 
Ikz 



-zxi/jqk 



-I- ichi^b+l--^j [V(n ^)x^pqk]^hu;qkz.x^pqk, 



cfikz 



rV'gk 



n{x,y) 

+ izchla- 



chAj 



(50) 



-V'qk 



1 



-2-,) (^(^ 



The operators applied to zx'i/jgk and i/jqk (in the left- 
hand side) show different longitudinal components (third 
term in each equation). At the same time, the gauge 
transformation (|46p implies transversality of spinor func- 
tions {zipqk—Q), such that these longitudinal terms cancel 
out by adjusting the gauge transformation parameters a 
and b : 



a={l + 2n)/2n, 6 = (1 - 2/x)/2^, 



(51) 



where /i is the magnetic permeability (a=§ and 6=—^ in 
the optical spectrum range). 

The vector cross product zxipqu in the first equation 
(150)) is the result of spin operator— Sz=i2X applied to 
the spinor rpq^. In the paraxial gauge transformation 
considered here (with Z'0=O), the spinor functions are 
invariant under the operator s^. Therefore, substituting 
parameters (1511) and taking a cross product of z and first 
equation, one obtains in ([50)) two identical eigenproblems 
that read 



n 



2mQ 

\n / \n[x, y) \n / 



(52) 



'lpqk = hu)ql^1pqk 



where mo—hkz/c{n~^) is the effective mass and pj^ = 
—iti\7± is the momentum operator in the lateral {xy 
plane) direction. The third term in the Hamiltonian 
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is the effective crystal potential induced by variations of 
the (complex) refractive index n{x,y). 

This is the stationary Schrodinger equation for photons 
in dissipative dielectric lattices of coupled microcavities 
or parallel waveguides. The orthogonality relationship 
can be established between the concomitant partners of 
biorthonormal set of its solutions [see Eq. ([50)1 ] 

(v^S'iV'S')- y'(tA;'k'-'^.k)dV=5(k^-ki)v,. (53) 

In the particular case of parallel dielectric waveguides, 
the conventional waveguiding theories have established 
the mode orthogonality relationship in terms of the elec- 
tric and magnetic field components. The biorthonormal 
set of solutions in mirror-patterned structures is shown 
to satisfy such relationship [Sec. |TT1 Eq. ([57|) ]. How- 
ever, the paraxial gauge transformation in periodic di- 
electric structures [Eq. differs from the one in mirror- 
patterned microcavities with uniform dielectric material 
in the cavity [Eq.(fT8|)]. Therefore, it is crucial to verify 
that the gauge transformation pS)) and the biorthonor- 
mal orthogonality relationship (j53p are in agreement with 
the well-established results of conventional waveguiding 
theories. 

An electromagnetic wave associated with the concomi- 
tant partner "i/jgk [in (|53p ] propagates in the same struc- 
ture as the wave (j46|) but with the complex conjugated 
refractive index and impedance: 



e,k=(^")*5(Z-'^)*-t/.,k, 
hgk=(^"'')*^(^'')*-2x^,k. 
Respectively, the inverse transformation reads 

^5k=(^°)*^"\^"")*-e,k, 
zx^,k=(^"'')*^"\^'')*-V. 



(54) 



(55) 



Substituting -zx [zx'^^-k.] [Eq.(l55l)] and Vgk [Eg. pT)) ] 
in the orthogonality relationship ([53]) , we obtain 



+^b J z.[zx(£-^e,k)]^^^^||^d^r^ (56) 

-^a I z.[(r'Vk')*xz]^^^^^d2r^. 

The second and third terms in the right-hand side of this 
equation are null. The first term coincides with expres- 
sion in Eq. ([M]) . yielding the orthogonality relationship 
[see Eq. (IMl)] 



z-[h*,k,xeqk]d^r_L= -6{k^-k'j_) 6, 



q'q- 



(57) 



The fields E^'k' and H^'k' of the wave associated with 
the concomitant partner '0g'k' are defined by relation- 



ship (|44| with the complex conjugated dielectric func- 
tion. Taking this fact into account, one can verify that 

(58) 

where integration runs over the entire structure of par- 
allel dielectric waveguides (or cavity-unfolded array of 
microcavities.) The expression (|58p corresponds to con- 
ventional orthogonality relationship between the modes 
of parallel dielectric waveguides. This result verifies the 
theoretical treatment presented in this paper. 

For low-order photonic bands, the structure parame- 
ters satisfying the paraxial approximation conditions can 
be defined more precisely. Thus, for a particular photonic 
band, the conditions ()48|) expressed in terms of intraband 
matrix elements read 



< 



I / \ I S 5 



(59) 



where the relationship Z{x,y)=fi/n{x,y) is taken into 
account. Being applied to the low-order bands, these 
expressions yield the relationships [see also the dis- 
cussion in Sec. Ill C| . 



IV. RESULTS AND DISCUSSION 

A. Generalized Hamiltonian 

Combining the results of Sec. [TTland lllll [Egs. (P5)l and 
([52)1 ] , we obtain the Hamiltonian for entire subclass of 2D 
structures characterized by paraxial light propagation in 
the direction normal to periodic lattice plane 



rriQC 



1 

n 

rriQC 



Pi. 

2mn 



i\ ch 

n/ 2L 
1 ' 



n{x,y) 



(60) 



■05k = ?it^gk'0gk- 



The photonic state wave function i/j^k (spinor) is related 
with the corresponding electromagnetic wave via paraxial 
gauge transformation 



qk 
.^9k , 



l+r;(z) 



<k\ 



(61) 



where operator £ was introduced in Eq. (jl9|) and param- 
eters a and h were obtained in Eq. (|5ip . 

The effective crystal potential is defined by the third 
and fourth terms in the left-hand side of Eq. ([60]) . In the 
case of 2D arrays of coupled microcavities, it takes into 
account the effects of mirror refiectivity patterning, as 
in metal-patterned VCSEL arrays^' (Fig. [T]) , and dielec- 
tric material variations, as in periodically etched VCSEL 
structures^i^ (Fig. [2]). Such simple expression for the 
effective potential was obtained by unfolding the cavities 
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along the optical axis [z-axis in Fig. [T]Jd)] and represent- 
ing the standing optical modes in the form of propagating 
Bloch waves (|6T|) in equivalent 3D structure. 

In a periodic array of microcavities, the longitudi- 
nal component of wave vector (kz) is fixed by the self- 
repetition condition at the cavity roundtrip. At the 
same time, Eq. (j60|) cannot reproduce the cavity reso- 
nance condition. Since Eq. (|60p takes into account a 
phase shift at the cavity mirrors, kz has to be evaluated 
from the roundtrip condition in a cavity with perfectly 
refiecting mirrors. In particular, kz—2'K/L in the case 
of one- wavelength microcavities {L=X/n). In addition, 
analyzing a symmetry of the group of k, the z-axis non- 
reciprocity of equivalent cavity-unfolded structure has 
to be taken into account (see Sec. HT])- Obviously, for 
an array of parallel dielectric waveguides or antiguides 
{e.g., photonic band gap fibers;^ photonic liquid crystal 
fibers^), these restrictions of the model do not apply. 
The effective potential is uniquely defined by the refrac- 
tive index variations [fourth term in the Hamiltonian of 
Eq.dSni)]- 

For most important (in practical applications) low- 
order bands, and within the range of parameters limited 
by conditions ([^T|) and the Hamiltonian is suitable 
for structures with high contrast of refractive index vari- 
ations {e.g., etched arrays of pillar microcavitieSf^ ho- 
ley photonic crystal fibers^). The class of 2D photonic 
crystal materials encompassing valid solutions of the non- 
Hermitian Hamiltonian eigenproblem (|60p is further illus- 
trated below with several structure examples. 

Arrays of semiconductor microcavities (n=3.5) defined 
by mirror refiectivity patterning (FiglT]) and operating at 
the optical wavelength A^l fim typically employ lattices 
of the pitch A^5 /im. These parameters assume that 
^2^10-3 [Eq. (gH)] and for any lattice fill factor FF, 
the Hamiltonian can be applied to the structures with 
mirror refiectivity contrast 6R up to 10~^. Note that in a 
typical VCSEL array with mirror reflectivity patterning, 
the refiectivity contrast is in the range from ~10~'' to 
-10-3. 

In periodic dielectric structures exhibiting high- 
contrast variations |(5n|>l (arrays of micropillars^* or 
holey photonic crystal fibersi^), the paraxial Hamilto- 
nian (|5D|) applies in the two opposite cases of lattice cell 
parameters [see Eq. dH])], at FF < X/2A\Sn\ (low fiU fac- 
tor) or at FF > l-X/2A\Sn\ (high fiU factor). Thus, for 
a 5 /im-pitch array of deeply etched semiconductor mi- 
crocavity pillars {n—3.5, \Sn\—2.5) operating at ~1 /im 
wavelength, low-order photonic bands can be treated by 
the model in the case of lattices with fill factor FF<0.2 
or FF>0.8. 

The first working silica photonic crystal fiber^" has 
a cladding material consisting of 300 nm air holes ar- 
ranged in a hexagonal lattice of 2.5 jim pitch. The re- 
fractive index of fused silica varies from 1.55 to 1.44 in 
the wavelength range of 0.2 — 1.5 nm. The lattice fill 
factor and the contrast are thus FF~0.02 and |(5n|'--^0.5. 
Due to the small size of air holes, the average refrac- 



tive index of the structure is close to that of fused silica 
[(l/n)-^ = {l+\Sn\)/{l+FF\6n\)r^l.5 in Eq. For 
such holey photonic crystal fibers, the Hamiltonian (pO]) 
is accurate throughout the entire optical transparency 
range of fused silica, from ultraviolet (A— 200 nm) to in- 
frared (A— 1.5 fim) wavelengths. 

In the case of arrays of parallel dielectric waveguides 
(or antiguides) with low-contrast refractive index vari- 
ations (|(5n|<Cl), the Hamiltonian (|60|) is valid for any 
fill factor of the lattice, in the optical wavelength range 
from A~4|(5n|A [Eq.dH])] to A--0.2nA. (The long wave- 
length range is limited by the paraxial approximation 
condition ^^1 with the critical value of £—0.1 .) The sil- 
ica photonic band gap fibers reported in Ref. ISg (n— 1.5, 
(5n=0.015 and A=6 fim) satisfy the paraxial approxima- 
tion conditions in the visible (A>360 nm) and near in- 
frared regions of the optical spectrum, up to the upper 
transparency edge of fused silica (A— 1.5 iim). 

B. Biorthonormal solutions in lattices with 
inversion symmetry 

The Hamiltonian in Eq. (j60p is independent of the spin 
variables. Therefore, all states are two-fold degenerate 
by spin (polarization), yielding thus degeneracy of the 
biorthonormal spinors ^pq■k^{t,r±) and xpq-ki{t,r±) asso- 
ciated with the eigenvalues ?iWgkT = ^^qkl ■ (Arrows indi- 
cate the spin direction. )^^ This fact allows the Eq. (|60p to 
be transformed into a scalar eigen problem with respect 
to the amplitudes of positive-spin (negative-spin) com- 
ponents of the spinors ^^)qk^_l{t,r±)=^^pqit^{t,r±)\'\ ,[) and 

■^9kT,i(i,rj_) = V'9k(i,rj.)|T,i) ■ 

In the case of 2D photonic lattices exhibiting inver- 
sion symmetry {e.g., square or triangular lattices), the 
biorthonormal system of lattice-periodic functions 

^,kTa(t,r±) = e-'"*+^''---u,kT,i(r±), (62) 
^,kTa(t,ri) = e-^"**+^''-'--u,kT,i(ri) (63) 

can be readily obtained by applying the PT transforma- 
tion (time reversal followed by coordinate inversion) to 
Eq. ([SO]) and noting the degeneracy ?ia;qkT ~ ^^gki ■ This 
has much in common with the Kramers degenerac y^^i^^ 
of a single-electron Hamiltonian in lattices with inversion 
symmetry. The invariance of the single-electron Hamil- 
tonian under the PT transformation yields degeneracy 
of eigenvalues i?gkT = ^-gki associated with the states of 
opposite spin. 

In the Hamiltonian (j60p . the degeneracy on spin vari- 
able {hLOqki;=huiqki) IS caused by the fact that the Hamil- 
tonian contains no spin operators. On the other hand, 
the Hamiltonian itself is not invariant under the PT 
transformation. In the lattices with inversion symme- 
try, its PT transform is uniquely defined by the result 
of the time reversal operation (complex conjugation fol- 
lowed by substitution —t t)r^^ which transforms H 
into H*=H+. 
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The PT transformation of the spinor functions (1621) 
reads 



Pr:i/;,kT.i(^,r^) = e 



—iLut-\-i'k.±r_i 



Uqk(l"±) 



„ — 'tcj* i+^k I r I 



(64) 



where transformation ahers the spin directio n^^i^^i^'^ (for 
s=l, It, i)*=— 1 1, T) in accordance with the phase con- 



vention of Refs. 43II49I ). Therefore, the PT transform of 



Eq. (|60p with the Hamihonian exhibiting the symmetry 
H*=H+ and H{—r±)—H{r±) yields, within the accu- 
racy of a phase factor, 



PT : iJt/>gkT,i(^.I"-L)=?it^gkT,i'09kTa(^.I"_L) 



where huj*^^-^—huj*^^. Comparison of with Eqs. ((27)) - ((28l) 
shows that the co-partners of biorthonormal set can be 
chosen as 



(66) 



'/'qkT,i(*'r_L)='0gk(i,r_L)|t,i> 

'^gkT,i(i'i"_L)=V'gkT,i("^'"'^-L)=V'gk(-i -r±)|T,i>- 

where Vgk(i, i'_l) and its PT transform ^qu{t,r±) = 
^p*^(—t,—r_l) are the scalar amplitudes of nonzero spinor 
component and its co-partner. In periodic lattices, these 
are the periodic Bloch waves with plane wave envelopes 
related via the PT-transform [see Eqs. 

In (j66p . the relative phase of co-partners lAgkt.i and 
■'/'qki-i is set by the orthogonality relationship PT|) . This 
relationship evidences that an arbitrary phase factor e*" 
can also be introduced in both equations , modifying 
thus the absolute phases of co-partners but preserving 
their relative phase shift. The fact that the co-partner 
phases have to be the same (at t=0) can be seen from the 
following observation as well: There should be no differ- 
ence between the co-partners in the case of Hcrmitian 
Hamiltonian. Therefore, the functions '0gkt.i(ij r_L) and 
'4^qk].lit,'r±) coincide in the limit lm(wqk)^0. 

Eq. (|66p defines the biorthonormal system of solu- 
tions in the case of lattices with inversion symme- 
try. It effectively transforms (|60p into a scalar Hamil- 
tonian eigenproblem with respect to the amplitude 
ipq]g_{t,rj_) of nonzero spinor component. The orthog- 
onality relationship for the biorthonormal set of scalar 
functions '0gk(i,r_L) and ipqk{t,r±)—ilj*y^{—t,—r±) follows 
from Eq. (PT|) : 



{■4>q'-k\i^qk) 



ipq'-k{-t,-rj_)il;qi,{t,r±)(frj_^Sq>, 



(67) 



Finally, note that in the case of significant difference 
between reflection coefficients of a mirror for s- and 
p-polarized waves in microcavities ( | lni?s/i?p|^^^ in 
([T])) or large polarization anisotropy in periodic dielec- 



tric lattices {\nx 



^ ), the spin-orbit coupling 



effects set in, rendering invalid the scalar approximation 
(|66p . Another example of spin degeneracy removal, for 
which the scalar approximation ((66|l is unsuitable, is the 
Coriolis-Zeeman splitting of photonic energy bands in 
nonpermanent gravitational fieldi^ 



Biorthonormal plane wave expansion in square 
lattices 



For square lattice structures (Sec. lIVPp . it is more 
convenient to express the biorthonormal set of spinor 
functions (1661) in the Cartesian coordinates basis, in the 



(65) 

form of functions xp= 



For a first rank tensor 



xp representing a state of spin s—1, the relationships — 
^(=^)=i(^i_i-V;i,_i)/V2 and V^''^ = (V'i,i+V'i,-i)/V2 pro- 
vide a transformation between the Cartesian coordinates 
and |s, rris) functions bases. (For paraxial photonic states 
considered here, the spinor component '(/''■^^=~*'0i,o is 
null and therefore not indicated explicitly.) 

The degenerate spinor functions (|66p can be repre- 
sented as linear combinations of positive- and negative- 
spin states 



, -iilpqk]~tpqkl) 7 -j('0gkT-'0gki) 
Wqk,it = 7= , Wqk.it — — 



tpqkS'- 



V2 

tpqk1:+1pqkl 
V2 ^ 



V'gk.j 



V2 

V2 



(68) 



yielding the biorthonormal set of solutions associated 
with the X- and y-polarized states of electromagnetic field 



'Jqk,<i~ 



rqki-t-r±) 



i^qki-t^-^l-) 



(69) 



Here, the spinor components are expressed in the 
basis of Cartesian coordinates. Like the functions 
(|66)) . degenerate spinor functions rpqk,x and i/'gk.y 
{^^qk.x—^^qk,y—^^qkT,i) couvcrt Eq. ([50)) into an eigen- 
problem with respect to scalar amplitudes ^pqk(t,r±), 
which assume the orthogonality relationship (|67p . 

For nonzero spinor components, the orthogonally po- 
larized states (IMl) show equal distributions |V'gk(i, r_L)P 
associated with the energy fiux along the z axis. In exper- 
iment, such states are observed as orthogonally polarized 
modes showing indistinguishable intensity patternsi^i^ 

The stationary Schrodinger equation (|60p is solved here 
using a biorthonormal plane wave expansion of the peri- 
odic crystal potential and Bloch functions ([55)) . In the 
stationary case, the time evolution of wave functions can 
be omitted, yielding the biorthonormal set of scalar am- 
plitudes V'ijk(rj_) and V'gk(rj_)=V'Jk(~''J-)- The relation- 
ship between co-partners of the set implies complex con- 
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jugated coefRcients of expansion: 

Mr±) = ^e^'^"'"" E C'9k(G) exp(zGr^), 

G 



These are the two complex conjugated matrix eigenprob- 
lem equations. Within the accuracy of complex eigen 
values and normalization condition for expansion coefR- 
cients, the eigenproblem ([50]) is thus converted into the 

(70) 

usual form encountering in conventional OPW expansion 
method. 



G 



where G is reciprocal lattice vector in the x-y plane. The 
difference in expansion coefficients in non-Hermitian and 
Hermitian Hamiltonian cases can be appreciated by ex- 
amining the Parseval theorem for coefRcients Cgk(G): 



(^,'k|^,k) = E C«'k(G)C,k(G) = V, 

G 

ft can be seen that J2g Cgk(G) = l, as opposed to the 
usual expression J^g |C'qk(G)|^ =1 (Parseval's theorem). 
Note that for the lattices with inversion symmetry dis- 
cussed here, the two sums are consistent in the limit 
Im(wqk)— >0 (Hermitian Hamiltonian case). In this limit, 

V'gk— *''0(jk and the phases of wave functions can be ad- 
justed to obtain real expansion coefRcients Cqk(G). 

Thus, in photonic structures exhibiting inversion sym- 
metry of the lattice, the biorthonormal plane wave expan- 
sion series differ from the usual OPW series by normal- 
ization condition for the amplitudes of spatial harmon- 
ics. The Schrodinger equation ([50]) can be then readily 
converted into a matrix equation for Cqk(G) coefRcients. 
Furthermore, the inversion symmetry of the lattice as- 
sumes a simple relationship between the the matrix ele- 
ments of operators H and . Thus, the periodic crys- 
tal potential and its adjoint operator are represented by 
series of lattice harmonics with complex conjugated co- 
efRcients 



-moc 



= E exp (- iGr_L ) , 

G 

C/(ri)+=;7(ri)*-E^G exp(-zGri) 



1 



n(r_L 



G 

where V-g=Vg due to the symmetry of crystal potential 
U{-rj_) = U{rj_). Substituting ^ and ^ into the 
Schrodinger equation (|60p and its PT-transform, mul- 
tiplying by exp(— i(kj^ + G')rj^) and integrating over a 
lattice cell in the x-y plane, one obtains the matrix equa- 
tions 



E 

G 



[( 



TTIqC 



?i'(ki- 



2mn 



-huj„ 



qk Og'G 



-Vg- 



E[K 

G 



G' 



?i'(ki- 



C9k(G)=0 
-G^2 



2™?; 



- — huj 



gk 



JG'G 



G G' 



C;k(G)=0. 



D. Band structure of square-lattice arrays of 
optical microcavities or parallel waveguides 



To study the light propagation behaviour in paraxial 
(71) photonic structures incorporating loss and gain distribu- 
tions, the non-Hermitian Hamiltonian (|60p is used here 
to analyze the structures with simple lattice symmetry 
and cell topology. Its numerical solutions are reported 
for square-symmetry lattices depicted in FigsH] and [21 
The results apply both to the structures deRned by mir- 
ror reflectivity patterning (e.i?., metal-patterned VCSEL 
arraysii) and to the ones exhibiting periodic variations of 
refractive index {e.g., etched VCSEL arrays^^ or matrices 
of antiguided VCSEL microcavities^). 

All matrices of microcavities considered here have sim- 
ilar lattice cell topology indicated in Figs[T]and[51 In the 
case of cavities with mirror reflectivity patterning, the 
position of optical microresonators is deRned by high- 
reflectivity square pixels separated by low-reflecting cav- 
ity mirror domains forming a grid pattern. The flU factor 
of such square lattice is the area ratio FF^a^ / , with 
a being the square pixel width and A being the lattice 
pitch. In similar dielectric lattices deflned by refractive 
index variations, the position of microcavities is set by 
square dielectric waveguide (antiguide) cores embedded 
in the background of the cladding material. Obviously, 
the same expression for the fill factor can be used to 
characterize such lattices of dielectric waveguides, with a 
being the waveguide (antiguide) core width. 

The Hamiltonian ([SO)) and the gauge transformation 

(72) 

(j6ip utilize equivalent, cavity-unfolded (3D) representa- 
tion of microcavities. As discussed in Sec. IIV Al the 
longitudinal component of propagation vector k is de- 
fined by condition kz—2'iT/L (in case of one-wavelength 
microcavities). Furthermore, by virtue of the dissipative 
effects rendering the opposite z-axis directions nonequiv- 
alent, the group of k in Eq. (I6ip contains only symmetry 
operations preserving the z-axis direction. 

In the particular case of square lattice, the cavity- 
unfolded 3D structure (2L-periodic in the z-axis direc- 
tion) has a tetragonal symmetry (A^2L) with symmor- 
phic space group Tq associated with £'4/1 point subgroup 
of rotations42i^i^ The reciprocal lattice is of the tetrag- 
onal symmetry as well with the first Brillouin zone (BZ) 
being of the rectangular prism shape [Fig[T](e)]. How- 
(73) ever, the crystal z-axis nonreciprocity implies that rota- 
tions and reflections of the D^h group altering the z- 
axis direction are not allowed. Therefore, for square 
arrays of microcavities, the group of k has a reduced 
symmetry characterized by C41, point group rotations 
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FIG. 3: Photon energy in square-lattice array of microcav- 
ities defined by mirror reflectivity patterning with the phase 
contrast of 51nii = ilO~^ (solid curves). Dashed curves in- 
dicate band structure of empty lattice (Ini? = 0). Other 
parameters are A = 5 /im, FF — 0.5, n = 3.5, L — 266nm 
(one- wavelength cavity optimized for A = 960 nm). 



in the A and T points of the BZ [at k=(0,0, /cz) and 
(±^,±j,fcz), respectively] and C2v symmetry in the Z 
point [at k= (±f,0,fc,) or (0,±^,fcz)]. 

FigsOlin] show the results of band structure computa- 
tions for square lattices of microcavities with the pitch 
A = 5 fim and lattice cell fill factor FF — 0.5. The 
photon energy dispersion and lifetime energy broadening 
curves along the high-symmetry lines A-Z-T-A in the BZ 
are obtained here using 25 lattice harmonics in Eg. ((73)) . 
In numerical simulations, the cavity structure is assumed 
to be optimized for 960 nm wavelength (operation wave- 
length range of GaAs/AlGaAs VCSEL structures with 
InGaAs quantum wells in the optical gain region). 

Figs. [3] and [J] detail the impact of mirror patterning 
contrast 6lnR — S\R\ + i6f and take contributions of 
the periodic phase (Sip) and amplitude variations 
of the mirror reflectivity into account. 

Fig. [3] shows the energy dispersion curves for low-oder 
photonic bands in the case of a lattice defined by phase 
modulation of mirror reflectivity with parameters 5|i?|=0 
and 6(p=10~'^ (black curves). The reflectivity at pixel po- 
sitions is i? = 1 while the periodic pattern i?(rj_) is intro- 
duced by a phase shift —S(p at reflections form the cavity 
mirror domains forming a grid. Comparison with the 
band structure of empty lattice (Fig|31 gray curves) indi- 
cates a blueshift of the optical modes due to contribution 
of phase-advancing domains of the grid. The main out- 
of-phase mode T5 (doubly degenerate by polarization)^^ 
has maxima of probability amplitude \ip\ located at the 
pixel positions and zeros located at the grid. Therefore, it 
shows the smallest energy shift. The contribution of the 
grid is particularly pronounced for Tg states exhibiting 
distribution with maxima located at the cross points 
of the grid and zeros at the pixels. 

Numerically calculated wave functions of these states 
(not shown in the figures) are in good quantitative agree- 
ment with the results^S obtained by means of group the- 
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FIG. 4: Photon energy (a) and lifetime energy broadening (b)- 
(d) in array of microcavities with refiectivity patterning con- 
trast SlnR oi 10~^ (black curves) 10~^ -|-ilO~"^ (gray curves) 
and 10~^ -I- zl0~^ (light gray curves). Other parameters of 
the structures are listed in the caption of Fig[3] 
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oretical analysis.' The intensity patterns of the opti- 
cal modes associated with these states are considered in 
Refs. [21II23 both theoretically and experimentally. 

Even for such low-contrast reflectivity patterns as the 
ones considered here, a complete 2D band gap can be 
opened in the lateral direction, between the T5 and Z2 
states. Thus, in Fig[31 the gap in the spectrum of optical 
modes is of 0.1 meV width. Since the effect of the grid on 
T5 states is small, the width of the energy gap is mostly 
defined by the blucshift of Z2 states of an empty lattice. 

Wave functions of the Z2 and states in the first BZ 
have zeros only along one crystalline direction (either x 
or y-axis directions) and differ by the position of proba- 
bility amplitude {ipl maxima. Along this lattice direction, 
the Z3 states are localized to the pixels while Z2 states 
are localized to the grid, like, respectively, and 
states. However, wave functions of the states originating 
from the Z points in the first BZ show large probability 
amplitude oscillations (and localization) along one lattice 
direction, as opposed to the states from the T point with 
wave functions oscillating in both lateral crystalline di- 
rections. For completeness of the discussion note that the 
wave functions of the lowest energy states in the A point 
(doubly degenerate by polarization Ai states) show no 
zeros of probability amplitude \ip\. Therefore, due to the 
effect of the grid pattern, the Z3 {Z2) states exhibit an 
intermediate blueshift as compared with the energy shifts 
of the Ai and T5 states (Ai and Tg states, respectively) 
of an empty lattice. 

Fig. [3] shows other states in the BZ (Z3 states) ex- 
hibiting small energy blueshift, which is comparable to 
that one of the states with wave functions localized 
to the pixel positions. The Z'^ states originate from the 
next nearest equivalent Z-points of reciprocal lattice at 
(±f,±2f,/c^) and (±2f,±f,fc^). Therefore, in addi- 
tion to the unidirectional oscillation features of the Z- 
states in the first BZ, their wave functions show large 
oscillations of probability amplitudes in the second crys- 
talline direction as well. The probability densities of the 
Z'2 states thus have a better overlap with the pixels of 
reflectivity pattern, which explains the smaller blueshift 
energy of Z'^ states as compared to their counterparts 
in the first BZ. In Fig. [31 the energy broadening of the 
bands due to the optical cavity loss has not been taken 
into account. 

Figd] illustrates the effects of optical loss distributions 
in the structures with the same lattice parameters as 
in Fig[3] (A = 5 /im and FF—0.5). All structures in 
Fig. [4] have the same amplitude reflectivity patterning 
|i?(r_L)| and differ only by the phase contrast of the pat- 
tern. The amplitude reflectivity |i?(rj^)| is 0.999 and 
0.989 for the pixels and grid domains of the pattern, re- 
spectively, yielding the amplitude contrast of (5|i?|=10~^. 
Note that the pixel reflectivity corresponds to the cavity 
loss of 0.2%, which is in the range of optical losses in a 
typical VCSEL structure. As in the case of the struc- 
tures in Fig. [21 in Fig. [H reflections at pixel domains of 
reflectivity pattern introduce no additional phase shift 



into the cavity roundtrip phase accrual of optical modes. 
Variations in the photon energy and lifetime broadening 
dispersion curves in Fig. [H are thus introduced by dif- 
ferent phase-advancing shifts at reflections from the grid 
domains in these structures. 

At no phase modulation of reflectivity pattern [6ip — 
0], the energy band structure is close to the one of an 
empty lattice but also shows a set of new peculiar features 
[Figs[ll(a), black curves]. Thus the degeneracy in energy 
of photonic states in high symmetry points of the BZ is 
partially removed (the states are doubly degenerate by 
polarization) and a set of partial flat bands appears in 
the Z- and T-points of the BZ. These features of energy 
dispersion curves are uniquely deflned by the dissipative 
effects in the structure. 

The energy broadening curves (or loss-dispersion 
curves) of optical modes in this structure are shown in 
Fig. m (b). (Only the losses related with the disperssive 
features of the patterned cavity mirror are taken into ac- 
count.) The doubly degenerate photonic states T5 have 
the lowest cavity loss (and lifetime energy broadening). 
They are of particular interest since they deflne the main 
lasing modes in coupled laser arrays i ^^i^^ The Z'^ states 
deflne the next lowest-loss modes in the high-symmetry 
points of reciprocal lattice. Finally, the states Tg are as- 
sociated with the highest-loss optical modes. Using the 
same considerations as in the case of energy dispersion 
curves in Fig. [H one can readily explain these features 
of the loss dispersion curves in terms of the overlap be- 
tween photonic state wave functions and high reflectivity 
pixels. 

Introduction of phase variations at the array grid does 
not affect the optical losses of modes associated with T5, 
Z'^ and states [see Figs[2 (b), (c) and (d)]. Fig[l 
shows the energy- and loss-dispcrsion curves calculated 
for the phase variations contrast Sip of 10~^ (gray curves) 
and lO^'^ (light gray curves). Note that in the last 
case d(p=S\R\. For small phase variations S(p<6\R\ (gray 
curves), the structure of energy dispersion bands is close 
to that one in the case of pure amplitude modulation of 
reflectivity pattern (see Fig[31(a), gray and black curves). 
However, for Sip^S\R\ (light gray curves), it approaches 
the energy band structure in the case of pure phase mod- 
ulation of mirror reflectivity (FiglSl black curves) and 
exhibits a band gap between T5 and Z2 states. 

A common feature of the loss-dispersion curves in 
Fig. [4] is the opened photonic band gap in the loss 
domain^ (or in the domain of photon lifetime in the cav- 
ity). Thus, it is impossible to excite an optical mode 
showing longer cavity lifetime than in the Tg state. By 
properly varying the lattice cell flU factor^Siii^ this ef- 
fect has been used to implement photonic crystal het- 
erostructure wells capable to conflnc photonic envelope 
wave functions to the regions of lower band gap mate- 
rial. The conflned states show usual features with cosine 
envelope functions in the well and exponential tails in 
the cladding material. The dissipative photonic crystal 
materials have been used so far in such photonic crystal 
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Fill Factor, FF 

FIG. 5: Double photonic crystal band gap: Band edge en- 
ergies of the Ts and Z2 states (left axis, black curves) and 
optical loss of the T5 state (right axis, gray curve) as a func- 
tion of the pattern fill factor FF in array of microcavities with 
reflectivity patterning contrast JlnJ? = (1 x 10"^. Other 
parameters of the structure are listed in the caption of Fig. [S] 
The double photonic crystal band gap is opened in the range 
of 0.16<FF<0.66. 

heterostructures show no band gaps in the photon energy 
domain. 

For a structure with complex parameter of reflectivity 
patterning contrast 8^p—6\R\—Vd~^ (the same parameters 
as in Fig. 01 light gray curves). Fig. [5] shows variations 
of the two photonic band gaps (in the photon energy 
domain and in the optical loss spectrum) as a function 
of the lattice cell fill factor FF . The lowest loss state 
T5 (gray curve, right axis) defines the band edge in the 
cavity loss domain. The gap below this edge is opened at 
any fill factor of the lattice. In the range of lattice cell fill 
factor 0.16—0.66, the energy of Z-2, state is higher than the 
band edge Ts, such that a second band gap originating 
at the band edge T5 exists in the photon energy domain, 
in parallel with the gap in the photon lifetime domain. 

The notion of double photonic crystal band gap il- 
lustrated in Fig. [5] opens new possibilities for tailoring 
photonic envelope wave function and controlling quan- 
tization features of confined photonic states in photonic 
crystal heterostructures. Thus, the eigenvalues of the 
Hamiltonian ((60)) at a photonic crystal heterostructure 
barrier assume the dispersion relationship 

fl k I 4 ?i k I r) 

nu,A^^r^ = 'hi^B^^r^. (74) 
Am A ■^nT'B 

where indexes A and B distinguish photonic crystal ma- 
terials at the heterostructure barrier and the parabolic 
band approximation is used for both photonic crystal 



materials. For complex band edge parameters TvjOa,b and 
effective masses m^.s, this condition assumes that both 
propagation constants k_L_^ and kj^^s are complex, in- 
dependent of particular photonic crystal heterostructure 
configuration. 

Thus, for an A^-dimensional photonic crystal het- 
erostructure well, Eq. (f74|) envisages a possibility of 2N- 
dimensional confinement of photonic envelope wave func- 
tions by introducing quantization of both real and imag- 
inary parts of propagation vector kj^ = k^ -I- zk" of con- 
fined photonic states. For such states, k" 7^ even in 
the region of lower band gap material (at the well core), 
allowing the confined states to be excited at the ener- 
gies within forbidden energy gaps of the heterostructure 
materials (well core and barrier materials). 

In the numerical examples of band structure computa- 
tions presented here, photonic structures utilizing reflec- 
tivity patterning for definition of periodic crystal lattice 
are considered. The effective crystal potential of these 
structures is governed by the third term in the Hamil- 
tonian of Eq. ((60|) and the magnitudes of its matrix ele- 
ments are bounded to the second-oder perturbations (of 
the relative order ~ ■C^). Therefore, in considered here 
case of paraxial light propagation with photon energy of 
1.3 eV and ^^~10^, the Hamiltonian (160]) yields accurate 
estimates of photonic bands splitting up to 1 meV . At 
the same time, for dielectric lattices defined by periodic 
variations of refractive index, the effective crystal poten- 
tial [fourth term in the Hamiltonian] is bounded to the 
first-order perturbations (~ C)- the case considered 
here {fiLo — 1.3 eF, 1^^0.03), the Hamiltonian ([60]) allows 
energy dispersion curves with band structure splitting up 
to 40 meV to be analyzed. 

V. CONCLUSION 

In this paper, a simple non-Hermitian Hamiltonian 
formalism is developed for a subclass of two dimen- 
sional photonic crystal structures characterized by parax- 
ial light propagation in the direction normal to the lattice 
plane, like in the case of coupled microcavity arrays or 
microstructured photonic crystal fibers. It allows the op- 
tical loss or gain distributions to be taken into account 
in the band structure analysis and envisages the effect 
of double photonic band gap opened both in the photon 
energy and lifetime domains. Predicted novel features of 
optical mode behaviour at double photonic crystal het- 
erostructure barriers with band edge discontinuities in 
the energy and lifetime domains offer new possibilities 
for photonic crystal applications in optoelectronic devices 
and integrated photonic circuits. 
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